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I  Introduction 


A  facility  providing  temporally  short  pulses  is  extremely  useful  for  investigating 
the  limitations  of  optical  detectors,  or  signal  processing  networks.  For  instance,  two 
picosecond  (ps)  pulses  have  a  500  GHz  bandwidth  and  can  therefore  be  used  in 
experimental  systems  designed  for  operation  one  to  two  orders  of  magnitude  faster 
than  those  presently  in  use.  This  report  describes  the  fiber-grating  pair  optical  pulse 
compression  set-up  at  the  Photonics  Center  in  Rome,  NY  which  compresses  110-120 
picosecond  pulses,  with  a  1.319  micron  {pm)  wavelength,  down  to  1-2  ps.  The 
exact  temporal  width  of  these  pulses  is  measured  by  an  autocorrelator. 

It  is  easy  to  compress  pulses  with  a  wavelength  <  1.32  (pm)  if  they  have  the 
following  qualities  a)  their  spectral  bandwidth  is  larger  than  the  inverse  of  their 
temporal  pulse  width  and  b)  the  instantaneous  frequency  varies  linearly  across  the 
pulse  (i.e.  it  has  a  linear  chirp).  In  a  fiber-grating  pair  compression  stage  the  light 
pulse  is  first  coupled  into  a  fiber.  Two  processes,  which  will  be  explained  briefly, 
transform  the  pulse  so  that  it  has  the  above  two  qualities.  Self-phase-modulation 
(SPM)  broadens  the  pulse  spectrally  satisfying  (a),  and  group- velocity  dispersion 
(GVD)  broadens  the  pulse  temporally.  When  the  optimum  fiber  length  for  the  given 
input  pulse  is  chosen,  the  output  pulse  has  a  linear  chirp.  At  this  point  the  pulse  is 
sent  through  a  dispersive  delay  line  with  anomalous  GVD,  which  exactly 
compensates  for  the  fiber  imposed  linear  chirp.  When  the  pulse  emerges  from  the 
delay  line,  which  in  our  case  is  a  grating  pair,  the  "red"  wavelengths  at  the  leading 
edge  of  the  pulse  are  delayed  just  enough  to  let  the  "blue"  wavelengths  at  the  trailing 
edge  of  the  pulse  catch  up.  The  net  result  is  a  large  temporal  compression  of  the 
pulse. 

The  report  begins  with  a  very  brief  theoretical  background  on  pulse  compression 
and  then  gives  the  step-by-step  procedure  for  setting  such  a  system  up.  Following 
the  pulse  compression  segment  is  a  short  description  of  the  autocorrelator  which  was 
built  to  measure  the  temporal  width  of  these  pulses. 

II.  Theoretical  Background 

The  next  several  sections  follow  closely  the  presentation  in  parts  of  Chapters  1, 
2,  3,  4,  and  6  of  Nonlinear  Fiber  Optics  by  G.  P.  Agrawal  (Academic  Press,  1989). 
First  the  basic  propagation  equation  of  a  pulse  in  a  fiber  is  set  forth.  This  equation 
contains  terms  responsible  for  both  temporal  pulse  broadening  (GVD),  and  spectral 
broadening  (SPM).  Additionally,  using  this  equation  it  is  easy  to  see  why  a  linear 


1. 


chirp  becomes  imposed  on  the  pulse. 


A.  The  Propagation  Equation 

The  electric  field  of  a  pulse  propagating  down  a  fiber  can  be  written 
E  (-,t)  =  ^{F(x,y)  A(z, t)  exp[i(/90z  -  w0t)]  +  c.c.}  (1) 


Here  F(x,y)  is  the  transverse  mode,  A(z,t)  is  the  slowly  varying  amplitude  of  the 
pulse  envelope  which  enters  the  fiber  at  z=0,  w0  and  are  the  pulse  envelope 

carrier  frequency  and  carrier  wavenumber,  respectively,  and  c.c.  stands  for  complex 
conjugate.  By  including  nonlinear  polarization  source  terms  in  the  wave  equation,  it 
can  be  shown  that  the  slowly  varying  amplitude  of  the  pulse  envelope  A(z,  t)  satisfies 


dA  a  dA 
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d*A 
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+  f  A  =  h\A\2A 


(2) 


The  coefficients  and  p2  take  into  account  the  dispersion  of  the  wavenumber  p(u) 
as  a  function  of  frequency,  that  is 


£(“>)  «  P0  +  (u  -  ui0)pl  +  |(w  -  w0)2^  +  ... 

Physically  p2  is  related  to  the  pulse  group  velocity  vg,  p2  =  — ,  and  P2  is  related  to 

the  GVD,  p2  =  Q  anc*  1  *n  ab°ve  equation  take  into  account 

absorption  losses  in  the  fiber,  and  nonlinearities  respectively  (a  =  n^^cA^)  where 
c  is  the  speed  of  light,  n2  is  the  intensity  dependent  part  of  the  refractive  index,  and 
Arff  is  the  effective  core  area). 

Equation  (2)  can  be  rewritten  in  normalized  units.  First  define  a  time  scale  in 
the  reference  frame  of  the  moving  pulse  which  is  normalized  to  the  initial  pulse 
width  T0, 

t  -  z/v,  X 
To  "  T0 


Also,  normalize  the  pulse  amplitude  to  the  incident  peak  power  P0 


') 


A(z,  r)  =  [P^  exp(-az/2)  u(z ,r) 

N 


The  propagation  equation  (Eq.  (2))  becomes 


.  3u  sgn(02)  32U 
ldz~  2Ln  dr* 


where  sgn(/?2)  =  ±1,  depending  on  the  sign  of  /92,  LD  =  j^-,  and  Lnl  =  — .  LD 

and  are  length  scales  over  which  dispersive  or  nonlinear  effects  become 
important  in  the  pulse  evolution.  In  what  follows  each  term  in  the  RHS  of  Eq.  (4) 
will  be  considered  separately  also,  since  we  are  interested  in  the  "normal"  dispersion 
regime  of  optical  fibers  we  use  sgn(/J2)  =  +1. 

B.  Group  Velocity  Dispersion  (GVD) 

The  fiber  can  be  considered  as  a  linear  dispersive  medium  by  setting  7  *  0  in 
Eq.  (4).  Then  Eq.  (4)  becomes 


•  _u  1  o  <92u 
dz  2P*  arr* 


The  Fourier  Transform  of  Eq.  (5)  is 


which  has  the  solution 


ti(z,w)  -  11(0,  w)  exp 
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So  at  a  point  z  in  the  fiber 
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Equation  (6)  shows  that  at  a  point  z  each  spectral  component  is  the  same  as  at  z=0, 
except  for  a  phase  factor  e‘*  where  <f>  =  /92w2z/2.  This  phase  factor  not  only  changes 
the  input  temporal  pulse  shape,  but  gives  it  a  frequency  chirp  also. 

Gaussian  Pulse  Example 

Consider  a  Gaussian  input  pulse  (z=0)  which  has  the  form 

u(0,T)  =  exp 

Using  Eq.  (6)  we  see  that  the  linear  dispersive  fiber  changes  its  shape  to 
U(Z’T)  =  T2  -  °i/$2z  CXP  ( 2(Tg  T  i^Z)  ] 

Notice  that  T2  —  T2  =  T2  |l  +  (recall  z/LD  -  |£2|z/T2).  Since  T2  >  T2  then 

the  pulse  height  is  reduced  and  the  pulse  becomes  temporally  broader.  Also,  by 
writing  u(z,T)  as  an  imaginary  number 

u(z,T)  =  ur(z,T)  +  UjjJz.T) 


=  |u(z,T)|  e^(*'T) 


we  find  that  the  instantaneous  frequency  across  the  pulse  varies  linearly  with  T 

Aw  =  -  =  - 2(-Z/Zo) - T 

3T  T2(l  +  (z/Ld)2) 

The  pulse  now  has  a  linear  chirp  as  a  result  of  the  GVD.  At  the  temporal  center  of 
the  pulse  (T=0)  the  instantaneous  frequency  is  w0,  but  at  a  leading  or  trailing  portion 
of  the  pulse  (T*0)  the  instantaneous  frequency  varies  from  w0  by  6w  as  given  in  the 


above  equation. 


C.  Self  Phase  Modulation 

Equation  (4)  shows  that  the  propagation  equation  also  has  a  term  proportional  to  the 
normalized  intensity,  |  u  | 2 .  To  see  the  effects  of  this  term,  we  momentarily  "turn 
off"  the  GVD  effect  by  setting  P2=0.  The  propagation  equation  becomes 

fr  «  ~ -  exp(-az)  |u|2  u  (7) 

az  Lnl 

which  has  the  solution 

u(z,T)  »  u(0,  T)  exp(i^,L(z,T))  (8) 

where 

^NL  (^T)  =  |u(0,T)|2  (9) 

Here  zrff  is  an  effective  fiber  length  that  takes  into  account  absorption  loss,  that  is 
zeff  *  ~a  0  -  exp(-atz)),  and  ^  z^  =  z. 

From  Eq.  (9)  we  see  that  the  instantaneous  optical  frequency  varies  across  the 
temporal  profile  of  the  pulse  from  its  central  value  w0.  The  amount  it  varies  by,  at 
T#0  is 

Aoj(T)  -  w(T)  -  w(T  -  0) 

<10) 

On  the  other  hand  Eq.  (8)  shows  that  as  the  pulse  propagates,  its  magnitude  does 
not  change;  i.e.  there  is  no  temporal  shape  change.  However,  the  phase  of  the  pulse 
does  vary  with  propagation  distance  z  and  intensity  |u|J,  as  is  shown  explicitly  in 
Eq.  (9).  The  consequence  of  this  self  phase  modulation  is  that  as  the  pulse 
propagates  the  spectral  component  at  each  point  in  the  pulse  varies  by  an 
increasingly  larger  amount  from  the  component  w0  at  the  center  of  the  pulse. 


5. 


Therefore,  under  the  right  conditions  an  incident  pulse  will  undergo  significant 
spectral  broadening  within  a  fiber. 

In  the  above  presentation,  either  GVD  or  SPM  was  considered  separately, 
however  both  effects  are  always  present  and  there  is  an  interplay  between  them. 
For  instance,  because  SPM  generates  a  broader  spectrum,  GVD  causes  the  pulse  to 
spread  out  temporally  more  than  if  there  were  no  SPM.  This  can  result,  in  a  pulse 
with  a  rectangular  temporal  profile,  and  linear  chirp.  Figure  3  shows  an  example  of 
this. 

The  relative  importance  of  GVD  versus  SPM  can  be  determined  by  rewriting 
the  propagation  equation  one  more  time  as 


•  2  3u  1 

it  d(z/z0)  2  dr 2 


N2  e~ai  ju  1  u 


where 
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■K 
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f-D  >  T  = 
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T 

*  o 


and 


N* 


Ld  _  kP0T2 

lnl  iSr 


GVD  dominates  for  N*  «  1  and  SPM  for  N2  >>  1.  At  1.319  /im  typical  values  of  7 

and  are  7  =  ■■■  ,  |/?2|  *  P5*.  For  a  100  MHz  modelocked  YAG  laser 

Wkm  km 

typical  values  of  P0  and  Tft  are  Pw  =  100  w,  T0  =  100  ps.  So  N2  =  10®  and  SPM 
dominates  as  would  be  expected  since  1.319  is  so  close  to  the  zero  dispersion 
wavelength  7S. 


/ 


D.  Grating  Pairs 

A  grating  pair  with  the  proper  configuration  can  act  as  a  delay  line  with 
anomalous  dispersion.  The  effective  GVD  parameter  is 


where 
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4jtc  b0 

3<:  *  wgA2cos2  8 o' 

Here  b0  is  the  center  to  center  spacing  of  the  grating  pair,  A  is  the  grating  period, 
and  all  other  quantities  are  either  as  defined  earlier,  or  defined  in  figure  1.  When 
the  proper  slant  length  b0  is  used,  the  grating  pair  compensates  for  the  fiber  induced 
linear  chirp  and  the  pulse  emerging  from  the  fiber-grating  pair  is  compressed. 

III.  Fiber-Grating  Pair  Compressor  Design 
A.  Useful  Formulas 

For  a  pulse  with  a  given  temporal  width  T0,  peak  power  P0  and  wavelength  A,  it 
is  useful  to  have  formulas  which  guide  in  the  fiber  compression  design.  The  needed 
quantities  are  a)  optimum  fiber  length  and  b)  grating  pair  separation.  It  is  also 
useful  to  know  c)  the  expected  compression  factor  T^T,.. 

a)  Optimum  fiber  length  zopt 
zopt  *  I.6/N  z0 

b)  Slant  length  b0  is  derived  from 

ac  K6 
T2  “  N 


so 


w2  d2  cos2  8ot 


4/rc 


UL  J2 
N  0 


c)  Compression  factor  -=? 


it-M 


or 


Te- 


L6  T 
N  0 
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If  fiber  lengths  z  are  used  which  are  much  shorter  than  zopt,  the  following 
formulas  a;,  used 

]3  fcng  d2  cos*  7  1  16  ^oTi 
0  84  (  4* c  J  N2  z  0 

j  „  _ L _ 

c  _  1  +  0.9N2  (z/z0) 


For  example,  typical  cw-modelocked  parameters  for  a  pulse  at  1.319  (im  are  T0  = 

1 10  ps,  P0  =  100  W.  Typical  parameters  for  a  dispersion  shifted  fiber  are  |02|  =  20 
ps2/km  Tf  =  20/Wkm.  So  z0  =  950  km,  N  =  1.1  x  10+3  and  using  the  above 
equations  for  z  ~  zopt,  we  get  zopt  =  1.4  km.  For  a  grating  with  1800  lines/mm  and 
using  Tf0  =  45°  b0  =  3  m. 


B.  Alignment  Procedure 

There  are  only  two  steps  to  aligning  a  fiber-grating  pair  pulse  compressor.  Step 
I  is  to  select  the  fiber  length,  and  step  2  is  to  set  the  slant  length  b0.  From  the 
preceding  example  zopt  =;  1.4  km  and  b0  =;  3  m.  These  are  approximate  values. 

Choosing  the  right  fiber  length  amounts  to  choosing  a  fiber  length  that  results  in 
an  output  pulse  with  a  rectangular  spectral  and  temporal  shape.  Therefore,  a 
spectrum  analyzer  and  sampling  scope  are  helpful,  but  not  absolutely  necessary.  As 
an  alternative  to  adjusting  the  fiber  length  from  its  approximate  value  to  optimum 
length  for  a  given  P0,  the  fiber  length  can  be  kept  constant  and  the  inDut  power  P0 
can  be  varied  either  by  using  a  A/2  plate  and  polarizer,  or  a  circular  neutral  density 
filter. 

Here  is  the  alignment  procedure. 

1.  Couple  into  and  out  of  the  fiber  with  a  20X  and  10X  microscope  objective, 
respectively.  Use  a  power  meter  at  the  output;  a  >  50%  throughput  is  adequate. 

2.  Optimize  the  power  through  the  fiber.  This  is  done  by  noting  at  what  power 
stimulated  Raman  scattering  (SRS)  just  begins  to  appear,  and  can  be  observed  in 
3  ways: 

a)  Bounce  the  output  pulse  off  a  grating  and  look  at  its  first  order  spectrum 
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using  an  IR  viewer.  The  Raman  line,  which  is  due  to  a  330  cm"1  phonon  mode  of 
silica,  will  appear  separate  from  the  band  centered  around  1.319  pm. 

b)  Look  at  the  output  spectrum  using  a  spectrum  analyzer.  The  Raman  line  is 
at  ~  1.405  fim. 

c)  Look  at  the  output  pulse  using  a  sampling  oscilloscope.  The  Raman  pulse 
will  lead  the  main  pulse. 

A  rule  of  thumb  is  that  ~  5%  of  the  total  intensity  should  be  in  the  Raman  line. 
Keeping  this  in  mind,  P0  can  be  increased  if  the  fiber  is  shortened. 

3.  Set  the  grating  spacing;  use  a  double  pass  configuration  so  that  b0  can  be 
charged  without  transversely  translating  the  grating  pair  output  path.  Adjust  b0 
for  the  minimum  pulse  width  by  monitoring  the  temporal  width  using  an 
autocorrelator.  It  may  be  necessary  to  coarsely  adjust  the  grating  spacing  using 
a  sampling  scope  and  fast  detector  for  the  initial  adjustment,  and  then  use  the 
autocorrelator  for  the  final  optimization. 

IV.  Autocorrelator 

The  autocorrelator  is  a  modified  Michaelson  interferometer.  A  time  variable 
delay  is  introduced  by  sending  both  beams  through  a  rotating  glass  block  ~  7  mm 
thick.  One  beam  is  displaced  vertically  by  a  right  angle  prism  before  the  two  beams 
are  recombined  in  the  second  harmonic  doubling  crystal.  The  flat  retroreflecting 
mirror  is  mounted  on  a  one-dimensional  translation  stage  that  is  useful  for 
calibrating  the  measurement  screen  of  the  oscilloscope  displaying  the  autocorrelator 
trace. 

V.  Experimental  Results 

Figures  2  and  3  show  actual  experimental  traces  of  the  temporal  and  spectral 
profile  of  a  1.319  pm  wavelength  at  different  stages  of  a  pulse  compression  setup. 
Figure  2  shows  the  initial  pulse,  which  has  a  temporal  width  of  101  ps  and  a  <lnm 
bandwidth.  At  the  output  of  the  fiber  (figure  3)  it  has  broadened  both  spectrally 
and  temporally  due  to  the  nonlinear  effects  discussed  earlier.  The  spectral  content  of 
the  pulse  has  broadened  to  3.5  nm  and  its  intensity,  as  a  function  of  wavelength,  has 
the  is  fairly  even  distribution  desired.  The  factor  of  two  temporal  broadening  is 
comparable,  but  not  equal  to  the  broadening  factor  predicted  by  the  above  equations. 
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This  is  not  surprising,  since  the  above  equations  art.  actually  intended  as  design 
guides  rather  than  stringent  design  rules.  More  careful  comparisons  of  theory  and 
experiment  must  be  done  using  beam  propagation  techniques  (see  Agrawal's  book  for 
an  overveiw  of  such  techniques).  Presumably  in  figure  3b  the  long  wavelength 
spectral  spectral  components  of  the  pulse  are  at  the  front  end  of  the  pulse,  the  short 
wavelength  components  are  at  the  back,  and  there  is  a  linear  variation  in  between. 
This  could  be  checked  by  an  elaborate  setup  such  as  a  streak  camera  measurement 
where  first  the  sampled  pulse  is  sent  through  a  dispersive  (nonchirping)  mechanism, 
such  as  a  spectrometer.  It  could  also  be  checked  by  taking  a  cross  correlation 
measurement  of  the  fiber  output  pulse  with  an  already  characterised  compressed 
pulse.  In  this  second  method  a  diode  array  at  the  output  of  a  spectrometer  would 
detect  the  different  (upconverted)  components.  The  second  method  is  doesn't 
require  a  streak  camera  so  is  much  easier,  but  it  is  also  of  course  redundant;  if  you 
have  a  well  characterised  compressed  pulse  why  characterise  for  the  sake  of 
compression  a  second  pulse  of  the  same  wavelength.  For  these  reasons  it  is 
preferred,  and  more  importantly,  adequate  to  rely  on  the  empirical  techniques 
mentioned  in  the  III  B,  to  determine  if  the  fiber  output  pulse  has  the  proper  spectral 
preparation. 

In  the  setup  which  produced  the  traces  of  figures  2  and  3  it  was  sufficient  to 
have  an  initial  grating  spacing  which  was  within  a  factor  of  3  of  the  final  optimum 
spacing.  This  resulted  in  about  a  20  ps  pulse  which  was  then  gradually  optimized  to 
about  3  ps  by  adjusting  slant  length.  Final  adjustment  to  the  observed  2  ps 
autocorrelation  pulse  was  obtained  by  <2  changes  in  the  pulse  launch  intensity.  This 
measured  pulse  corresponds  to  an  actual  pulse  width  of  about  1.3  ps  once  the  1.5 
deconvolution  factor  (for  gaussian  pulses)  is  included. 

It  is  worth  mentioning  that  when  the  glass  block  is  symmetrically  situated  to 
generate  time  delays  in  both  beams,  as  in  this  case  (see  figure  4),  the  generated  time 
delay  versus  angle  change  is  linear  over  a  surprisingly  large  range.  This  is  confirmed 
by  a  simple  ray  diagram  calculation.  The  useful  range  of  the  autocorrelator  of 
course  depends  on  the  thickness  of  the  glass  block  used  as  a  delay  line.  Increasing 
the  autocorrelator  range  of  course  means  that  the  block  will  have  to  spin  faster  for 
live  monitoring  purposes.  This  puts  higher  demands  cn  the  mechanical  components 
the  autocorrelator,  and  increases  the  need  for  proper  dynamic  balancing.  A  smaller 
range  (thinner  block)  is  useful  only  once  the  system  is  near  optimization.  For  a 
system  producing  1-2  ps  pulses  a  useful  (i.e.  linear)  range  of  20  ps  seems  to  be 
adequate  once  the  above  tradeoffs  are  considered. 
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Appendix  A:  Useful  Formulas 


1.  Factorized  form  of  electric  field  in  a  fiber 


E(-’t)  =  i  {p(x’y)  A(z’tJ  exPt‘0V  -  “oO]  +  c.c.} 


2.  Propagation  Eq.  for  slowly  varying  part  of  E(r,  t) 


dA  .  a  dA  -a  d*A  .  aA  _  :~|a|2A 
-to  +01  ~to+1^  w+  2A  *  inr,A| 


a  is  absorption,  7  = 


n0w0 


cA 


eff 


where  n2  =  3.2  •  lO-20  m2/w 
w0  =  carrier  frequency 
c  =  speed  of  light 
A«fr  *  effective  fiber  core  area 
Also 

01  =  —  is  vg  =  group  velocity 

v« 

and 


0  . - j-  is  GVD 

v£  dw 


note:  pt  is  related  to  the  group  delay  dispersion  D(A) 


D  -  -  also  |D(A)|  -  D(A)  •  c  •  A 
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Figure  1.  Pulse  compression  setup 
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Figure  2.  Input  pulse:  101  ps 
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Figure  3.  Output  pulse  after  fiber: 

a)  Frequency  spectrum,  AX  =  3.5  nm; 

b)  Temporal  profile.  At  *  225  ps 


Parts  Lists 


The  major  components  of  both  the  pulse  compressor  and 
some  vendors,  are  listed  below. 


A.  Pulse  Compressor 

part 

vendor 

(2)  X/2  plate  @l.319um 

CV1 

polarizer 

CVI 

Fiber  coupler 

Newport 

Fiber  coupler 

Newport 

Microscope  objective 

Newport 

fi 

Newport 

Rotation  stage 

Control  Optic 

Diffraction  Gratings 

Spectron 

1.5  km  Dispersion  shifted  fiber  (1.55um) 

(3)  Mirror  mounts 

Newport 

B.  Autocorrelator 

part 

vendor 

Color  Filter 

Oriel 

f  =  50mm  lens 

Melles  Griot 

Beamsplitter 

Ealing 

Li  103  SHG  crystal  (in  cell) 

Inrad 

Mirror 

Melles  Griot 

autocorrelator,  along  with 


n 

F-1015LD-DJ 
F  916  T 
FL10B 
FL20 
RS37A 

2"  sq„  1200  gr/mm, 
1.319um  blaze 
Corning 
MM2 

ti 

51312 

01  LPX  108 
356139 

10  deg.  convergence 
angle 

02  MFC*  000.  A!  flat 
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Turning  prism 
(I)  1"  Mirror  mount 
(1)  2"  Mirror  mount 

(4)  iris  diaphragm  Edmund  Scientific 

Phototube  Hamamatsu 

Phototube  Housing 

15V  DC  power  supply  for  phototube 
2"  x  1"  x  3cm  BK-7  glass  flat 
Electric  motor  for  glass  block 
5V  Power  supply  for  electric  motor 


D30.263 
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MISSION 

°f 

Rome  Air  Development  Center 


RADC  plans  and  executes  research,  development,  test  and 
selected  acquisition  programs  in  support  of  Command,  Control, 
Communications  and  Intelligence  (C*I)  activities.  Technical  and 
engineering  support  within  areas  of  competence  is  provided  to 
ESD  Program  Offices  (POs)  and  other  ESD  elements  to 
perform  effective  acquisition  of  C9I  systems.  The  areas  of 
technical  competence  include  communications,  command  and 
control,  battle  management  information  processing,  surveillance 
sensors,  intelligence  data  collection  and  handling,  solid  state 
sciences,  electromagnetics,  and  propagation,  and  electronic 
reliability /maintainability  and  compatibility. 


